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gold and silver nanospheres illuminated by a tightly focused beam. Using a multipole-expansion 
approach we obtain an analytical solution to the scattering problem and provide insight on the 



{S) effects of focusing on the optical response. Large differences with respect to Mie theory are found 



especially when the nanoparticle supports quadrupole or higher-order resonances. 
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I. INTRODUCTION 

The optical properties of metal nanoparticles (NPs) [H [2] represent a topic of active 
research in several areas, like nano-optics [31 SI IS] , solid-state physics E], biology 0, 
[TUl [TT] and photonics [121 E3]- What makes them particularly attractive is their ability to 
sustain an electromagnetic resonance yet being much smaller than the incident wavelength. 
This peculiar feature originates from the existence of an electron-density mode that can 
couple to external radiation giving rise to a localized surface-plasmon-polariton resonance 
(SPR) [H, 12]. In the experimental practice, metal NPs have been often studied close to 
interfaces and under both wide-field and tight illumination [UJ [THJ US! EZ]- Nevertheless, 
since the NP is very small compared to the wavelength, the measured plasmon spectrum is 
usually only compared to theory based on plane wave (PW) illumination [1] . While the effect 
of supporting films and substrates has been extensively investigated [HI [191 EQJ, [2H [22] , to 
the best of our knowledge, there is no detailed study on the spectrum of metal NPs under 
tight illuminations [23, [21] ■ 

The electromagnetic problem of a linearly-polarized PW scattered by a spherical particle 
dates back to the work of Mie [25]. After the invention of the laser, the availability of narrow 
band and collimated beams has posed new theoretical issues on the scattering of light by 
particles. Morita et al. [26], and Tsai and Pogorzelski [2Z] were the first to study the case 
of spherical particles illuminated by a Gaussian beam whose waist is smaller or larger than 
the particle size. A decade later, Barton et al. [2H] extended the study from the far field to 
the near and internal fields. Other authors have also considered the scattering properties of 
objects under Gaussian illumination [221 EH EH E2]- In summary, all works have found that 
there is a significant deviation from the case of PW illumination if the particle size is larger 
than the beam waist. 

Less attention has been given to the case of a high-numerical-aperture (high-NA) beam 
incident upon a metal NP [221 121] • Torok et al. [23] arid Challener et al. [23] have described 
the high-NA beam by PW expansion [33J . Subsequently, the scattering problem was solved 
using classical Mie theory [23] or the finite-difference time-domain method [24J. Although 
the PW expansion provides an accurate description of the electromagnetic field at the focus, 
it does not provide an easily accessible physical insight on the effect of strong focusing. 
Moreover, in the mentioned studies, only selected wavelengths were considered and the 



modification of the plasmon spectrum was not taken into account. 

The aim of this work is to study the interaction of a high-NA beam with metal NPs in 
more detail and compare it with the case of PW illumination. By solving the scattering 
problem using a multipole-expansion approach [31], we gain intuition on the NP's optical 
response both in the far and near field. We show that even when the NP is smaller than 
the focal spot, the difference in the multipole content of a high-NA beam compared to a 
PW causes a modification of the plasmon spectrum, especially when the metal NP exhibits 
resonances of order higher than the dipolar one. 

The paper is organized as follows. In Sec. |TT] we present the analytical solution of gener- 
alized Mie theory [35] for a high-NA beam incident on a spherical NP placed at the focus. 
We also derive formulae for the scattering and extinction cross sections and for the average 



intensity enhancement in the near field. In Sec. |IH| we present and discuss results for gold 
and silver NPs with a 100 nm diameter. 

II. THEORY 

A. Generalized Mie Theory for a High-NA Beam 

The generalized Mie theory of a high-NA beam incident on a spherical NP placed at 
the focus is carried out by expanding the electromagnetic field in multipoles. The inset to 
Fig. [I] sketches the layout of the problem. A PW polarized along x and propagating along 
—z is focused by a high-NA aplanatic system [33]. It can be shown that the incident field 
components at the reference sphere in the image space are 



El, r = 0, (1) 

E lfi = E e' lkf a(6) cos 0, (2) 

E- h4> = -E e~ lkf a(6) sin 0, (3) 

in the spherical coordinates (r, 8, <fi). Here E Q is the PW amplitude and k is the wavevector. 
The parameter a(8) represents the illumination angular weighting factor and / is the lens 
focal length. The factor exp(-ikf) is added to ensure phase agreement with the electric 
field computed by PW expansion [33]. For a system satisfying the sine condition one finds 
a (9) = (cosfl) 1 / 2 [33]. 



The multipole-expansion coefficients of the incident field in the image space are obtained 
by enforcing the boundary conditions set by Eqs. (flip]) [34]. The <p dependence in these 
equations already tells us which multipoles can properly represent the incident field, namely 

Ei(r, 9, <P) = J2 [4N ell (r, 6, 0) + B,M oli (r, 9, <P)] , (4) 

i 

where the sum over I extends from 1 to infinity and N e ij, M D iz are the vector spherical 
harmonics defined in Ref. [U The explicit dependence of the electromagnetic field on (r, 9, <f>) 
is omitted in the following expressions for the sake of brevity. Furthermore, we choose the 
spherical Hankel function of the second kind (h\ (p)) [36] in the multipoles so that the 
incident field is an incoming wave. Using the asymptotic expression h\ (p) ~ iexp[— i(p — 
ttI/2)]/p and matching the far fields of Eq. Q to Eqs. ({l]|3]), we obtain the conditions 
B\ = —iA\ and 
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where P/(cos#) are associated Legendre functions J3U] and p = kf. Using the procedure 
described in Ref. |3H the expansion coefficients are found to be 
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The parameter a is the angular semi-aperture of the lens and is determined by the formula 
NA=rib sin a, where n^ is the refractive index of the background medium. The regular- 
ity of the electromagnetic field at the focus requires that Eq. Q contains also outgoing 
multipoles represented by spherical Hankel function of the first kind (h\ (p)) [36], with 
the same amplitude of the incoming ones, i.e. no sources at the origin [31]. That implies 
K (p) + ty (p) = 2ji(p), where ji(p) is the spherical Bessel function of the first kind [36J. 
We use the notation N^j and M^{ for multipoles with j;(p) and include the factor 2 into 
Ai to eliminate it from Eq. Q. 

When a spherical NP of radius a is placed at the focus, we have to solve for a boundary 
value problem involving the incident field Ei, the scattered field E s and the internal field 
E int . The sum of incident and scattered fields builds up the total field E tot = E ; + E s in 
the background region [T]. The symmetry of the problem implies that the NP will not excite 
multipoles except those present in the incident field. Therefore, we can formally write the 
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FIG. 1: Relative strength of the multipole coefficients for the expansion of a PW \E[/ E\\ and a 
high-NA beam {a = 60° and a = 30°) |A//Ai|. Inset: a PW polarized along x and propagating 
along —z is focused by a high-NA lens and scatters on a metal NP positioned at the focus, a is 
the angular semi- aperture of the lens. 
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where the notation N^j and M^ 1Z means ^ (p) in the multipoles since the scattered field 
must be an outgoing wave. The magnetic field H is simply obtained from Maxwell's equa- 
tions and from the properties of spherical vector harmonics jl]. 

Fulfilling the boundary conditions at the NP surface r = a leads to an expression for the 
expansion coefficients that reads 

mipi{mx)ip[{x) — ipi(x)ijj' l (mx) 



a, 



b, 
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mxi{x)4>[{rnx) — ipi(mx)x'i(x) ' 

where ipi(p) = pji(p) and Xi(p) — ph-i (p) are Riccati-Bessel functions, x = ka is the size 
parameter and m = k\jk is the relative refractive index pQ, if k\ is the wavevector inside 



the NP. Similar expressions hold for the coefficients of the internal field 

mxi(x)i)\(x) - mjji(x)x'i(x) 
■mxi{x)il)[{mx) — ipi(mx)x'i(x) ' 
d = mipi(x)x'i(x) - mxi(x)ipi(x) ,_ 

mipi(mx)x'i(x) — Xi{x)^[{Tnx) 



Notice that the coefficients of Eqs. (10 13) have exactly the same expression as the Mie 
coefficients for an incident PW pp. Once they are known, we simply use Eqs. ([7]{9]) to 
compute the electromagnetic field inside and outside the NP analytically 

B. Far- field Cross Sections 

The interaction of light with material particles is usually investigated by considering far- 
field cross sections. The scattering and absorption cross sections, C s and C a , represent how 
strongly a particle scatters or absorbs light, respectively, for a given incident intensity. The 
total amount of power that is removed from the incident beam via scattering and absorption 
is proportional to the extinction cross section C e = C s + C a p] . 

Since the radial component of every multipole vanishes in the far field, the scattered 
power W s reads 

W s = lim - / / Re {ErfHlt - E^H^} r 2 sin 0d0d0, (14) 

where E s = (E Str ,E St Q,E s>l p) and H s = (H Btr ,H Bt g,H s ^) represent the scattered field. After 
performing the integration over the angles and exploiting the orthogonality of the angle- 



dependent functions T\ = dPj 1 (cos 9) / d9 and iii = Pf (cos 9) /sin 9 [I], Eq. (14) simplifies 
to 

I 
where Z is the background medium impedance. Likewise, the extinguished power W e is 

defined as 

W e = lim - / / Re {E^H: e - E- lfi H* - E sfi H* + E B ,+H? e } r 2 sin#d0d0, (16) 
r ^°° % Jo Jo 

where E ; = (E- 1>n E^g, E^) and Hi = (H i<r , H ie , H^) refer to the incident field. After 
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Notice that Eqs. (15) and (17) have the same dependence on Mie coefficients as for a PW [T]. 
Since the intensity of a PW is homogeneous, the cross sections are easily obtained by dividing 
the scattered or extinguished powers by the illumination intensity. However, for a high-NA 
optical system the intensity at the focus is inhomogeneous, and introducing a cross section 
is not straightforward. Nevertheless, in order to be able to compare our results with that of 



a PW PQ, we define the cross section by dividing Eq. (15) and (17) by the incident intensity 

L '=^S^' ym -' is - (18) 

Here we compute the Poynting vector along —z and take its average over the NP xy circular 
section of area A = na 2 at the focal plane. 



C. Average Intensity Enhancement 

Far-field quantities like the cross sections do not provide information on the near field. 
To further compare the optical response of a NP illuminated by a high-NA beam with a 
PW, we would like to consider a parameter similar to the cross section as an estimate of 
the complex structure of the near field. Messinger et al. [37] have defined the near-field 
quantity Qnf by computing the scattering cross section at r = a, including the radial field 
component, and dividing it by the NP area. They also considered a similar quantity, called 
Qr, by taking only the radial component of the scattered field. They found that Qr is 
very close to Qnf because in the near-field the radial component is quite large. However, to 
study the near-field more precisely, one has to take into account both the incident and the 
scattered fields. Thus we evaluate the total field intensity |E tot | 2 , average it over the NP 
surface and then normalize it by the intensity at the focus to obtain an enhancement factor, 

K= . , w * ,2 / F r|Etot(a, M )| 2 sin^d0. (19) 

An |Ej(U,U,U)| Jo Jo 

The same average intensity enhancement can be computed for PW illumination to compare 
the two situations. It turns out to be more instructive if one looks at the radial and tangential 
components of the total field separately. Therefore, we split K into a radial K t and tangential 
K t average intensity enhancement. By taking advantage of the orthogonality properties of 



multipoles, these quantities have the analytical expression 

I 2 O/3/l i i\3 
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^ t = 16(^2. |^-im- [|ai1 |X;I + N lxi1 +w +m ~ (21) 

2Rfi{o»Xj}^-2Re{6,xj}^], 

where ^z(p) and Xi(/°) are calculated at the surface of the NP, that is p = ka. The coefficient 
A\ at the denominator derives from the normalization |E;(0,0,0)| 2 = 4|v4i| 2 /9. It can be 
shown that the same quantities for a NP illuminated by a PW have exactly the same form, 
except that Ai is replaced by E h where E\ = E i l (21 + 1)/ (1(1 + 1)) is the multipole-expansion 
coefficient for a PWpQ. 

When the NP is very small compared to the wavelength, only the first multipole con- 
tributes to the intensity enhancement both for a high-NA beam and a PW. By using the 
asymptotic expansions of Mie coefficients and Riccati-Bessel functions for p — > [lj, we 
obtain the following simple expressions, 

K T = - [4|a| 2 + 4Re{a} + l], (22) 

K t = ?[| a | 2 -2Re{a} + l], (23) 

where a is a quantity proportional to the NP polarizability per unit volume and is given by 

m 2 — 1 , . 

a = — . 24 

m 2 + 2 v ; 



Equations (22 23) contain contributions from the scattered field, the incident field and their 
interference. Similar formulae have been proposed to explain the electromagnetic contribu- 
tion to surface-enhanced Raman scattering in the limit of very small NPs [3H] • 

III. RESULTS AND DISCUSSION 

To understand the interaction of a high-NA beam with metal NPs, we choose gold and 
silver NPs with a diameter of 100 nm embedded in low-index lossless dielectric media, such 
as glass (n-b = 1.52) and air (rib = 1)- The NP size is assumed to be smaller compared to the 
focal spot, but sufficiently large to show deviations from PW illumination. Moreover, while 
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100 nm gold NPs exhibit mainly a dipolar SPR, 100 nm silver NPs can support higher-order 
SPRs. This allows us to investigate the role of SPRs in determining the optical response 
as a function of focusing. The analytical formulae of Sec. [TT] have been computed using 
Mathematica [39]. Since the size parameter for 100 nm NPs illuminated in the wavelength 
range from 300 nm to 900 nm is close to one, we include up to 15 multipoles in the field 
expansions. Including more multipoles did not bring a noticeable change in our results. For 
modeling gold and silver NPs we use the optical constants from Ref. HO) 

As it can be seen in Fig. [TJ if the beam gets focused more tightly, the dependence on 
higher-order multipoles becomes weaker compared to the case of PW illumination. Since 
the same multipole coefficients are present also in the scattered field, we expect that the 
excitation of higher-order SPRs will be strongly suppressed by high-NA beams. 

A. Far-field Cross Sections 

Figure [2] shows the scattering and extinction cross sections of a 100 nm gold NP in glass 
illuminated by a PW and a high-NA beam (a = 60°). Since the NP can support almost only 
a dipolar SPR, the difference between the two illumination conditions is small. The SPR 
peak is larger for a high-NA beam because the average incident intensity /; is smaller than 
that of PW illumination, as shown in the inset to Fig. [2^,. If we had chosen the intensity at 



the origin for the normalization of Eqs. (15) and (17), we would have not seen an increase in 
the dipolar peak. For smaller gold NPs, the difference becomes negligible since the average 
incident intensity is computed over a smaller area. As shown in the inset, for a given NP 
and background medium the average incident intensity depends on the wavelength and on a. 
The inset to Fig. [2b shows that the incident Poynting vector is indeed completely along — z 



across the NP, making the definition of average intensity according to Eq. ( 18 ) a reasonable 
choice. 

Figure [3] shows the scattering and extinction cross sections of a 100 nm silver NP in 
glass illuminated by a PW, a tightly (a = 60°) and a moderately (a = 30°) focused beam. 
Because the silver NP can support higher-order SPRs, the effect of focusing results in a 
strong suppression of the quadrupole mode at A ~ 420 nm and of the next higher-order 
mode at A ~ 390 nm. The latter is visible only in the extinction cross section because it 
mostly contributes to absorption rather than scattering. As in the case of gold, the dipolar 
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FIG. 2: Scattering (a) and extinction (b) cross sections of a 100 nm gold NP in glass illuminated 
by a PW and a high-NA beam (a = 60°). The inset to (a) shows the incident intensity averaged 
on the NP according to Eq. ( |18| ) and normalized with respect to the intensity at the origin as a 
function of the angular semi-aperture a for A = 450 nm and A = 650 nm. The intensiy at the 
origin is equal to that assumed for PW illumination. The inset to (b) plots the Poynting vector 
for the high-NA beam (a = 60°) at A = 450 nm with the 100 nm gold NP to scale. 



resonance at A ~ 530 nm is enhanced due to the average incident intensity. The difference is 
slightly larger than for gold because the resonance occurs at a shorter wavelength, where the 
focal spot and the average incident intensity are smaller (see inset to Fig. |2k). Therefore, 
by adjusting the focus one can manipulate the excitation of SPRs in NPs and the optical 
response can be significantly changed even if the NP size is smaller than the focal spot. For 
smaller particles the differences with respect to PW illumination vanish. 
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FIG. 3: Scattering (a) and extinction (b) cross sections of a 100 nm silver NP in glass illuminated 
by a PW and by a high-NA beam (a = 60° and a = 30°). 

B. Average Intensity Enhancement 

We now focus our attention on the near field to see how the field enhancement depends on 
the incident illumination. Figure [4] compares the average intensity enhancement of a 100 nm 
gold NP in glass. Since the gold NP exhibits mainly a dipolar SPR, also in the near field the 
difference between focused beam and PW is small. In contrary to Fig. |2j the dipole peak is 
not enhanced by a focused beam because the normalization is with respect to the field at 
the origin. In fact, a closer look at the data shows that the average intensity enhancement 
is larger for a PW than for a focused beam at the dipolar SPR. A small difference can be 
noticed also around A ~ 550 nm, where a weak quadrupole SPR exists. K r is significantly 
larger than K t as already found in Ref. |371 where only the scattered field was considered. It 



is interesting to notice that, as predicted by the terms in Eqs. (20) and (21), the radial and 
the tangential fields are maximal at different wavelengths. 

Figure [5] shows that for a 100 nm silver NP the average intensity enhancement does depend 
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FIG. 4: Average radial (solid lines) and tangential (dashed lines) intensity enhancements evaluated 
at the surface of a 100 nm gold NP in glass illuminated by a PW and a high-NA beam (a = 60°). 

on the strength of focusing. Like the far-field cross sections, the higher-order SPRs are weakly 
excited by a high-NA beam. Therefore, the enhancement is much smaller compared to PW 
illumination in the spectral region where the NP exhibits a quadrupole or a higher-order 
SPR. On the other hand, the enhancement associated with the dipolar SPR is left almost 
unchanged. Also for the silver NP, K v is larger than K t even though for higher-order modes 
the difference is less pronounced. It is also worth mentioning that while for the dipolar SPR 
K T and K t are maximal at different wavelengths, they reach their maxima at almost the 
same place for the higher-order SPRs. 

The radial enhancement for the dipolar SPR is red-shifted with respect to the dipolar 
peak in the far-field cross sections for both gold and silver [37]. The fact that in going 
from air to glass the SPRs responsible for the near-field enhancement are red-shifted is a 
well established result pp. However, modes higher than the quadrupole are not excited by 
the high-NA beam (see Fig. [5]). Indeed, Fig. [I] confirms that multipoles of order I larger 



than 2 are almost negligible in the beam expansion. According to Eqs. (23) and (22), in 



the asymptotic limit of very small NPs there is no difference between PW and focused 
illumination because both contain the electric dipole multipole. 
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FIG. 5: Average radial (solid lines) and tangential (dashed lines) intensity enhancements evaluated 
at the surface of a 100 nm silver NP in glass and air illuminated by a high-NA beam (a = 60°) (a) 
and a PW (b). 

C. Intensity Enhancement 



K r and K t provide us with a lot of information on the near field of metal NPs under 
different illumination conditions. However, since they are average quantities, they eliminate 
information on the spatial distribution of the enhancement. For this reason, we have to 
study the electric field pattern in the vicinity of the NP. We focus our attention on gold at 
the peak wavelength of the dipolar SPR obtained from K r , since the higher-order SPRs in 
the silver NP would be weakly excited by the high-NA beam. Moreover, we discuss only the 
radial component because the tangential parts do not exhibit a large enhancement. Again, 
it is worth mentioning that this quantity is normalized with respect to the intensity at the 
focus. 

Figure M depicts the radial intensity enhancement for a high-NA beam (panels (a) and 
(c)) and a PW (panels (b) and (d)) incident on a 100 nm gold NP in glass. As in the case of 
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FIG. 6: Radial intensity enhancement (log 10 scale) for a 100 nm gold NP in glass illuminated 
(along —z) by a high-NA beam (a = 60°) at A = 645 nm ((a) and (c)) and by a PW at A = 642 nm 
((b) and (d)). Cross cuts: xy plane ((a) and (b)), xz plane ((c) and (d)). The field is not computed 
inside the NP. 

K T we consider the radial component of the total field. Notice that the wavelength chosen 
for the high-NA beam and the PW are slightly different because the dipolar peaks in Fig. [4] 
are not exactly at the same position. We have checked that choosing the same wavelengths 
does not change the following results since the SPRs are quite broad. For the cross cuts 
corresponding to the focal plane (panels (a) and (b)) the field distribution looks very similar 
close to the NP surface, while a deviation occurs away from it because the incident field 
starts to be as strong as the scattered one. The enhancement is maximal along the x axis 
for both high-NA beam and PW. Also for the cross cuts in the xz plane (panels (c) and 
(d)) the patterns are quite similar in proximity of the NP, but differences occur at larger 
distances. This can be explained by the different field profile of a high-NA beam and a PW. 
As in the case of the far-field cross sections and the average intensity enhancement, there is 
no remarkable difference between a high-NA beam and a PW at the dipolar SPR. 

To have a closer look at the radial intensity enhancement of Fig. |6j we plot it at the NP 
surface as a function of 9 and for <fi = 0°, that is in the xz plane. Figure [7] shows both the 
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FIG. 7: Radial intensity enhancement as a function of 6 for = 0° at the surface of a 100 nm gold 
NP in glass illuminated by a high-NA beam (a = 60°) at A = 645 nm and a PW at A = 642 nm. 
The contribution only from the scattered field is represented by dashed lines. 

enhancement for the total intensity (solid lines) and for only the scattered intensity (dashed 
lines). First, it is seen that the constructive interference between incident and scattered field 
has a non-negligible contribution to the intensity enhancement. Second, the maximum for a 
high-NA beam and a PW are not at the same location. For a high-NA beam the peak value 
is reached at 9 ~ 100°, while for a PW it is reached at 9 ~ 105°. Moreover, the overall profile 
of the curves is slightly different. The enhancement for a high-NA beam is more symmetric 
with respect to 9 = 90° and broader than that for a PW. On the other hand, the maximum 
enhancement is larger for a PW than for a high-NA beam. We have checked that when the 
NP gets smaller these differences disappear and the maximum radial intensity enhancement 
occurs at 9 = 90° for both illuminations. 

IV. CONCLUSIONS 

Using a multipole-expansion approach we have studied the interaction of a high-NA beam 
with noble metal NPs placed at the focus and have compared the results to the case of PW 
illumination. The most important and general effect of focusing resides in the suppression 
of higher-order multipoles in the incident field expansion. Consequently, even if the NP 
supports high-order SPRs, they cannot be efficiently excited by a high-NA beam, as we 
have seen for the silver NP both in the far and near field. Therefore, care must be taken 
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when interpreting the plasmon spectra measured under high-NA illumination. 

In this work we have considered high-NA beams generated by linearly-polarized light. It 
could be interesting to extend these investigations to radially-polarized beams, which achieve 
a tighter focus jUJ S2] . These results might find application in areas like NP spectroscopy [HH 
El], optical data storage |l5j HE] and optical forces |Hj. 
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